A radiomics signature derived from CT imaging to predict MSI status and immunotherapy outcomes in gastric cancer: a multi-cohort study

Background Accurate microsatellite instability (MSI) testing is essential for identifying gastric cancer (GC) patients eligible for immunotherapy. We aimed to develop and validate a CT-based radiomics signature to predict MSI and immunotherapy outcomes in GC. Methods This retrospective multicohort study included a total of 457 GC patients from two independent medical centers in China and The Cancer Imaging Archive (TCIA) databases. The primary cohort (n = 201, center 1, 2017–2022), was used for signature development via Least Absolute Shrinkage and Selection Operator (LASSO) and logistic regression analysis. Two independent immunotherapy cohorts, one from center 1 (n = 184, 2018–2021) and another from center 2 (n = 43, 2020–2021), were utilized to assess the signature’s association with immunotherapy response and survival. Diagnostic efficiency was evaluated using the area under the receiver operating characteristic curve (AUC), and survival outcomes were analyzed via the Kaplan-Meier method. The TCIA cohort (n = 29) was included to evaluate the immune infiltration landscape of the radiomics signature subgroups using both CT images and mRNA sequencing data. Results Nine radiomics features were identified for signature development, exhibiting excellent discriminative performance in both the training (AUC: 0.851, 95%CI: 0.782, 0.919) and validation cohorts (AUC: 0.816, 95%CI: 0.706, 0.926). The radscore, calculated using the signature, demonstrated strong predictive abilities for objective response in immunotherapy cohorts (AUC: 0.734, 95%CI: 0.662, 0.806; AUC: 0.724, 95%CI: 0.572, 0.877). Additionally, the radscore showed a significant association with PFS and OS, with GC patients with a low radscore experiencing a significant survival benefit from immunotherapy. Immune infiltration analysis revealed significantly higher levels of CD8 + T cells, activated CD4 + B cells, and TNFRSF18 expression in the low radscore group, while the high radscore group exhibited higher levels of T cells regulatory and HHLA2 expression. Conclusion This study developed a robust radiomics signature with the potential to serve as a non-invasive biomarker for GC’s MSI status and immunotherapy response, demonstrating notable links to post-immunotherapy PFS and OS. Additionally, distinct immune profiles were observed between low and high radscore groups, highlighting their potential clinical implications. Supplementary Information The online version contains supplementary material available at 10.1186/s12885-024-12174-0.


Background
Gastric cancer (GC) ranks as the fifth most prevalent cancer and the fourth leading cause of cancer-related death globally, contributing to approximately 770,000 annual fatalities [1].While surgical resection followed by adjuvant chemotherapy is established as an effective treatment for early-stage GC, a significant proportion of patients are diagnosed beyond this stage due to the disease's insidious onset and rapid progression [2,3].Palliative chemotherapy serves as the mainstay treatment for unresectable or metastatic GC patients, yet its clinical benefits are very limited, resulting in a median survival time of less than one year [4,5].
In recent years, tumor immunotherapy has made remarkable strides in both research and clinical practice, offering a promising treatment avenue for patients with unresectable or metastatic gastric cancer [4,6,7].The combination of PD-1/PD-L1 inhibitors and chemotherapy has now become the standard first-line treatment for GC in major guidelines [8][9][10].Unfortunately, due to the strong heterogeneity of gastric cancer, only a minority of GC patients respond to current immunotherapy treatment [11].Thus, there is an urgent need to identify reliable biomarkers for screening GC patients likely to benefit from immunotherapy.
Microsatellite instability (MSI) stands out as the first pan-cancer immune biomarker approved by the Food and Drug Administration (FDA), with solid tumors exhibiting MSI-high (MSI-H) recommended for immunotherapy [12,13].The MSI-H phenotype generates numerous immunogenic neoantigens detected by the immune system, rendering MSI status a valuable clinical biomarker for checkpoint immunotherapy [14].Presently, MSI status assessment primarily relies on immunohistochemistry or polymerase chain reaction (PCR) analysis of specimens obtained through endoscopic biopsy or surgical resection [15].However, information on MSI expression status obtained postoperatively has limited impact on treatment planning prior to surgery.Moreover, the limited samples obtained through biopsy may not comprehensively reflect tumor heterogeneity, leading to falsenegative results (2.1-5.9%)[16,17].Additionally, biopsies and surgeries are invasive, time-consuming, expensive, and carry risks of complications, making repeated monitoring inconvenient [17].Hence, there is an urgent need to develop a non-invasive, reliable, and cost-effective method to identify MSI status.
Radiomics is a rapidly advancing field that utilizes advanced computational techniques to transform medical images, such as CT and MRI, into quantitative features, enabling the development of a signature for cancer diagnosis and treatment [18,19].While several studies have demonstrated the potential and significance of radiomics in evaluating the MSI status of gastric cancer, these investigations have not thoroughly delved into the clinical value of their radiomics models within cohorts of patients undergoing immunotherapy [20][21][22].Specifically, there is a lack of comprehensive validation regarding the efficacy of radiomics models in predicting immunotherapy outcomes.Moreover, the aforementioned studies did not delve into the potential biological value of radiomics models from the perspectives of the immune microenvironment and transcriptome.
Therefore, this study aims to fill the gaps in the current research landscape by establishing a non-invasive radiomics biomarker.This biomarker will not only identify MSI status, but will also undergo further validation in cohorts of patients undergoing immunotherapy, which will offer a comprehensive understanding of its practical feasibility.Moreover, we will analyze the immune microenvironment features of patients stratified into different radscore groups based on transcriptomic data, providing a deeper understanding of the predictive mechanisms of the radiomics model.

Patients
The study received approval from the local Ethics Committee (2022-KY-1447-002) in accordance with the Declaration of Helsinki.Written informed consent was waived by the Human Scientific Ethics Committee of the First Affiliated Hospital of Zhengzhou University due to the retrospective design of the study.
Patient selection process is depicted in Fig. 1.To develop the radiomics signature associated with MSI status, the study retrospectively enrolled 201 consecutive gastric cancer (GC) patients from Center 1 (the First Affiliated Hospital of Zhengzhou University, 2017-2022).MSI status was determined using PCR detection.The primary cohort was randomly divided into a training cohort (n = 142) and a validation cohort (n = 59) at a 7:3 ratio.In addition, two independent immunotherapy cohorts were employed to investigate the association of the radiomics signature with immunotherapy response and validate its prognostic value.The first cohort, ZZU cohort, comprised 184 GC patients treated at Center 1 between 2018 and 2021.The second cohort, SDU cohort, included 43 GC patients treated between 2020 and 2021 at Center 2 (the Second Hospital of Shandong University).Furthermore, one GC cohort (n = 29) from The Cancer Genome Atlas (TCGA) and The Cancer Imaging Archive (TCIA) databases was included in the study to further evaluate the immune infiltration landscape across different subgroups of the radiomics signature.Supplementary A outlines the inclusion and exclusion criteria for all cohorts, and Supplementary B provides information on the immunotherapy regimens used in the immunotherapy cohorts.

Detection of MSI status
MSI status was determined by performing PCR testing on tumor tissue obtained from surgical resection or biopsy.Five Bethesda microsatellite markers (BAT25, BAT26, D2S123, D5S346, and D17S250) were utilized to assess MSI status.Tumors with instability at two or more of these five markers were classified as MSI-H, while those without instability or with instability at only one marker were classified as MSS or MSI-L.In this study, we grouped tumors with MSI-L and MSS together as the MSS group due to their similar clinical significance.

CT acquisition and radiologic evaluation
CT scan parameters are detailed in Supplementary C. A consensus review of all images was conducted by two radiologists, LL and LP, who evaluated the clinical T stage and clinical N stage of each patient based on the AJCC version 8 TNM staging system.LL had 5 years of experience in abdominal CT imaging, while LP had 10 years of experience.Any discrepancies in radiologic evaluation were resolved through discussion and consultation with a third senior radiologist, GJ, who had over 20 years of experience in abdominal CT imaging.The radiologists were blinded to clinical information during the evaluation process.

ROI segmentation and feature extraction
Portal venous phase CT images of all patients were analyzed using the open-source software 3D Slicer (version 4.13.0;https://www.slicer.org/).Manual ROI segmentations were performed by outlining the tumor lesion's margin on the largest slice.To ensure intra-observer reliability, ZP performed the ROI segmentation and feature extraction process twice within a week for CT images of 50 randomly selected patients.For interobserver reliability assessment, two radiologists, ZP and YS, independently conducted ROI segmentation and feature extraction on medical images of the same 50 patients.ZP completed the remaining ROI segmentation independently.
To minimize the influence of different scanning schemes or equipment on quantitative radiomics features, the ROI segmented images were resampled to a voxel size of 1 × 1 × 1 mm 3 and discretized using a bin width of 25 for gray level discretization [19].A total of 851 radiomics features were extracted from each ROI segmented image.Following feature extraction, all radiomics features underwent standardized using the z-score for subsequent analysis.

Feature selection, model construction and evaluation
All 851 radiomics features were analyzed further based on the training cohort.To reduce the dimensionality of radiomics features and identify robust predictors of MSI status, we utilized a meticulous three-step procedure.Firstly, we calculated the inter-observer and intra-observer correlation coefficients (ICCs) for each feature and retained those with ICCs greater than 0.8 for further analysis.Secondly, a multivariate ranking method called maximum relevance minimum redundancy (mRMR) algorithm was applied to eliminate redundant and irrelevant features on the basis of a heuristic scoring criterion, and only the top ranked 200 features were retained.Lastly, the least absolute shrinkage and selection operator (LASSO) logistic regression algorithm was performed to choose the most valuable subset of features from the top ranked 200 features.The regular parameter (λ) of LASSO regression was chosen when the average mean square error was minimal by 10-fold cross-validation.Moreover, the most valuable subset of features was utilized to create the radiomics signature in the training cohort.
Clinical risk factors associated with MSI status were identified through univariate analysis with a significance threshold of p < 0.05.Based on the clinical risk factors and radiomics signature, we developed three models through multivariable logistic regression: the clinical model, the radiomics model, and the combined model.Discrimination performance was assessed by calculating the area under the receiver operating characteristic curve (AUC), and the results were compared using the non-parametric Delong test.The predictive accuracy of the radiomics signature was further evaluated in both the training and validation cohorts using calibration curves and the Hosmer-Lemeshow test.Furthermore, a nested 5-fold cross-validation approach was conducted on the selected radiomics features using logistic regression to further assess the performance of the radiomics signature.In the outer 5-fold cross-validation loop, we divided the data into training (4-folds) and test (1-fold) datasets to evaluate the performance of the models with an untouched test set.In each training step of the outer fold, an inner 10-fold cross-validation was applied using all the selected radiomics features to tune the hyperparameters and select the optimal logistic regression model.The nested 5-fold cross-validation process was repeated 5 times to ensure the robustness of our evaluation results.

Association between radiomics signature and immunotherapy outcomes
Patients' immunotherapy responses were categorized into four groups according to the RECIST 1.1 guidelines: complete response (CR), partial response (PR), stable disease (SD), or progressive disease (PD).The immunotherapy cohort patients were then grouped into two categories based on their treatment response: the CR/ PR group and the PD/SD group.The radiomics signature calculated a radscore for each patient in the immunotherapy cohorts, and its ability to predict a CR or PR response to immunotherapy was evaluated using AUC value.Additionally, the prognostic significance of the radiomics signature was evaluated through Kaplan-Meier survival curves and log-rank tests.Progression-free survival (PFS) was defined as the time from the initiation of immunotherapy to disease progression, while overall survival (OS) was defined as the duration from the initiation of immunotherapy to disease-related death or the last follow-up date.

Immune infiltration assessment
In the TCIA cohort, patients were categorized into a high radscore group (≥ the median) and a low radscore group (< the median).Utilizing mRNA sequencing data, the CIBERSORT algorithm estimated the abundance of 22 immune cell in the tumor immune micro-environment (TME).Furthermore, we assessed differences between the two groups concerning MSI status and immuneregulating factors, including co-stimulators, co-inhibitors, and other relevant factors.Immune regulators are known to play a crucial role in modulating the function of immune cells, exerting either anti-tumor or pro-tumor effects.

Statistical analysis
Statistical analyses were conducted using R software (version 4.2.2, https://www.r-project.org/).Continuous variables, if normally distributed, were reported as mean ± standard deviation; otherwise, they were expressed as median (lower quartile, upper quartile).Fisher's exact tests were applied to categorical variables, while appropriate tests such as Student's t-test or Mann-Whitney U test were used for continuous variables.ICC analysis, mRMR algorithm, and LASSO regression analysis were conducted using the "irr", "mRMRe", and "glmnet" packages, respectively.Receiver operating characteristic (ROC) curves were generated using the "pROC" and "ggplot2" packages.The nested 5-fold cross-validation was performed by applying "nestedcv" package.The prognostic value of the radiomics signature was assessed using Kaplan-Meier survival curves and log-rank tests.The "survminer" package was used to plot the survival curves, and the "IOBR" package [23] was conducted for immune infiltration analysis.The MSI status in the TCIA cohort was obtained using the "cBioPortalData" package, and a threshold of 0.4 was applied to distinguish between MSI and MSS [24,25].Statistical significance was defined as a p-value of less than 0.05.

Patient characteristics
A total of 457 GC patients were enrolled in this study.Demographic and clinical characteristics of the training cohort (n = 142) and validation cohort (n = 59) are summarized in Table 1.Within the patient population, 149 were male, 77 were aged 65 years or older, and the MSI-H positivity rate was 22.4%.
Clinical characteristics of the two immunotherapy cohorts are presented in Table 2.In the ZZU cohort, 132 patients were included, with 71.7% being male and 38.0% aged 65 years or older.Among these patients, 64 (34.8%) achieved CR or PR with immunotherapy, with median PFS and OS of 5.8 (2.9, 10.7) months and 10.8 (6.0, 16.2) months, respectively.The SDU cohort comprised 43 patients, with 81.4% being male, and none were aged 65 years or older.Among these patients, 18 (41.9%)achieved CR or PR with immunotherapy, and the median PFS and OS were 5.8 (3.7, 7.9) months and 9.5 (7.2, 14.3) months, respectively.
Baseline data for the TCIA cohort (n = 29) is provided in Supplementary TableE1.This cohort consisted of 29 patients, with 85.7% being male and 78.6% aged 60 years or older.The median radscore in this cohort was 1.45 (0.87, 1.84), and 14 (48.3%)patients were classified into the low radscore group.

Development and validation of the radiomics signature
Significant associations between tumor location and MSI-H expression were identified in both the training (p = 0.007) and validation cohorts (p = 0.015) through univariate analysis of clinical factors.To enhance robustness and eliminate redundancy, ICC analysis and the mRMR algorithm were employed to eliminate radiomics features.Subsequently, a 10-fold cross-validation LASSO algorithm (Fig. 2) was applied to develop a radiomics signature, comprising nine features.Supplementary Table E2 provides the coefficients of each feature in the radiomics signature.Clinical and radiomics models were independently created based on tumor location and radiomics signature.The combined model was developed by integrating tumor location and radiomics signature using multivariate logistic regression (Table 3).Detailed explanations for the three logistic regression models are presented in Supplementary Results.
Table 4 presents the diagnostic performance indicators, including accuracy, sensitivity, and specificity of the three models in both the training and validation cohorts.3e and f (p < 0.001).Finally, the results of the nested 5-fold cross-validation (Supplementary Results) indicated that the radiomics signature's performance remained reliable.

Clinical prognostic validation of radiomics signature in immunotherapy
The ZZU cohort exhibited a significant lower radscore (median:  4a-d.Furthermore, the results from the ZZU and SDU cohorts revealed that patients in the low radscore group had a significantly longer median PFS of 8.3 months (4.6, 12.8) and 7.4 months (4.9, 8.0) and a longer median OS of 14.0 months (10.4,18.7) and 12.2 months (9.0, 21.2), respectively, compared to those in the high radscore group.The high radscore group had a median PFS of 4.6 The results of Kaplan-Meier analysis revealed a significant association between the radiomics signature and both PFS (HR: 1.59 (95% CI: 1.12, 2.27), p = 0.009) and OS (HR: 2.00 (95% CI: 1.34, 2.98), p < 0.001) in the ZZU cohort (Fig. 4e and f ).Similarly, in the SDU cohort (Fig. 4g and h), the radiomics signature was significantly associated with both PFS (HR: 3.12 (95% CI: 1.31, 7.40), p = 0.003) and OS (HR: 2.51 (95% CI: 0.86, 7.35), p = 0.025).Importantly, patients in the low radscore group, predicted by the radiomics signature, experienced a significantly longer median PFS and OS compared to those in the high radscore group.These findings indicate that the radiomics signature can serve as a useful tool in predicting the survival benefit of immunotherapy in GC.

Immune infiltration in the high and low radscore groups
The immune infiltration landscape of both groups is depicted in Fig. 5.The results indicated an activated immune microenvironment in the low radscore group but a potentially immunosuppressive state in the high radscore group.Specifically, the low radscore group demonstrated significantly higher CD8 + T cell levels compared to the high radscore group (Fig. 6a).Although memory resting CD4 + B cells were lower, activated CD4 + B cells were significantly higher in the low versus high radscore group (Fig. 6a).Moreover, regulatory T cells (Tregs) were increased in the high radscore group compared to the low group (Fig. 6a).Additionally, the proportion of microsatellite instability (MSI) was notably higher in the low radscore group (Fig. 6b), along with increased TNFRSF18 expression (Fig. 6c).Meanwhile, HHLA2 expression was significantly elevated in the high versus low radscore group (Fig. 6d).

Discussion
The detection of MSI status is crucial to guide immunotherapy for GC patients, with MSI-H as a predictive biomarker [12].We developed and validated a non-invasive radiomics signature with good performance to predict MSI-H in GC.Importantly, our investigation revealed an association between the signature and immunotherapy response and outcomes in multi-center cohorts.Additionally, findings exhibited an activated immune microenvironment in the low radscore group, while the high radscore group showed immunosuppression.This underscores the biological significance of our signature in predicting immunotherapy efficacy.The prevalence of MSI-H in GC has been reported to vary from 5.6 to 33.3% in previous studies, possibly due to differences in detection methods.In this study, we used PCR as the gold standard for detecting MSI-H and observed an incidence rate of 22.4% in our study population.Compared to the deficient DNA mismatch repair (MMR) immunohistochemistry (IHC) method commonly used in most studies, PCR can more accurately reflect the MSI status.Furthermore, the IHC method is highly dependent on sample quality, and samples may gradually lose nucleic acids or proteins over time, leading to decreased detection accuracy [26].
Previous studies have demonstrated that certain clinical characteristics, such as age, gender, and tumor location, are significantly correlated with GC MSI [27].A recent meta-analysis [28] found that women have a significantly higher probability of exhibiting MSI than men (OR 1.57, 95% CI: 1.31 to 1.89; p < 0.001), and that MSI is significantly associated with those aged 65 years and older (OR 1.58, 95% CI: 1.13 to 2.20; p < 0.001) and upper GC location (OR 0.38, 95% CI: 0.32 to 0.44; p < 0.001).In this study, MSI-H expression was notably correlated with tumor location, consistent with previous research.The pathogenesis of GC varies depending on its location, leading to distinct expressions of phenotype markers, biological behaviors, and gene expression profiles [29,30].This may be the biological basis for the higher expression of MSI-H in gastric antrum cancer.However, age and gender did not demonstrate significant correlations with MSI, possibly due to a limited sample size and data bias.Moreover, MSI exhibited notable correlations with Lauren classification (intestinal subtype), TNM stage, and other pathological features [31].However, as the aim of this study was to establish a pre-treatment non-invasive prediction model, these factors were not taken into consideration.
In addition to clinical and pathological features, several studies have indicated that CT characteristics may hold value as predictors of MSI status.A retrospective study of 77 patients revealed that dMMR GC typically displays a lower stomach location, smaller tumor thickness and lymph node diameter, and fewer lymph nodes on CT imaging [32].Another study conducted by Wu et al., which analyzed clinical and pathological information from 114 patients with colorectal cancer, showed that multi-parameter analysis derived from single-source dual-energy CT can relatively accurately distinguish between MSI and MSS in colorectal cancer [33].Nevertheless, further large-scale studies are necessary to validate these findings.
Radiomics holds significant promise in predicting the MSI status of tumors [34][35][36].While the precise validation cohorts, respectively [22].However, they did not further confirm its clinical value in immunotherapy effectiveness and prognosis, nor did they deeply investigate its biological mechanisms.
In this study, we utilized the LASSO algorithm to construct a radiomics signature associated with MSI-H expression in GC.The top-ranking features, determined by their coefficient weights, include original_glcm_MCC, original_gldm_Large Dependence Low Gray Level Emphasis, and wavelet.LHL_glszm_Large Area High Gray Level Emphasis.These features offer insights into image texture complexity (MCC), gray level distribution (Large Dependence Low Gray Level Emphasis), and gray heterogeneity (Large Area High Gray Level Emphasis).The findings indicated that MSI-H gastric cancer potentially exhibits a more intricate texture structure, decreased gray heterogeneity, and a more uniform gray level distribution.Furthermore, our results demonstrated that the signature displayed promising performance, with an AUC of 0.851 (95% CI: 0.782, 0.919) in the training cohort and 0.816 (95% CI: 0.706, 0.926) in the validation cohort.These AUC values surpassed those reported in previous studies, underscoring the potential clinical utility of our signature.
Moreover, this study extended its investigation to assess the radiomics signature's predictive capacity for immunotherapy efficacy and prognosis.Notably, patients with a low radscore demonstrated a higher likelihood of achieving CR/PR following immunotherapy.The radiomics signature exhibited predictive AUC values of 0.734 (95% CI: 0.662, 0.806) and 0.724 (95% CI: 0.572, 0.877) for CR/ PR in the ZZU and SDU cohorts, respectively.Importantly, patients identified as low radscore through the radiomics signature displayed notably extended median progression-free survival (PFS) and overall survival (OS) compared to those in the high radscore group, suggesting its potential role in anticipating the survival benefits of immunotherapy in GC.Significant parallels can be drawn from Huang et al. 's extensive multi-cohort GC study [37], which established a CT-based radiomics score (RS) using 2272 patients.This study explored the correlation between the radiomics biomarker and the neutrophil-to-lymphocyte ratio (NLR) within the tumor immune microenvironment, including its link to prognosis and immunotherapy response in advanced GC. Huang et al. found that patients with lower RS (60.9% and 42.9%) exhibited substantially higher objective responses to anti-PD-1 immunotherapy compared to those in the higher RS group (8.1% and 14.3%).Moreover, Sun et al. integrated various solid tumor patient cohorts and developed an independently validated radiomics biomarker for tumour-infiltrating CD8 cells [38].They demonstrated its correlation with the tumor immune phenotype and its predictive potential for clinical outcomes in patients undergoing immunotherapy.These studies collectively underscore the micro-level role of radiomics features in predicting immunotherapy outcomes through the discrimination of the tumor immune microenvironment.
Considering this, the study conducted further analysis in the TCIA cohort to investigate the MSI expression status and immune infiltration in the two subgroups.The results consistently demonstrated a significantly higher proportion of MSI in the low radscore group, providing additional evidence for the strong correlation between the radiomics signature and MSI status.Notably, the immune infiltration analysis revealed compelling findings.The low radscore group exhibited significantly elevated levels of CD8 + T cells and CD4 + B cells activated compared to the high radscore group, indicating enhanced activation and functionality of these effector immune cells.On the other hand, the high radscore group showed notably higher levels of CD4 + B cells memory resting and Tregs compared to the low radscore group, suggesting a more immunosuppressive microenvironment with increased regulatory T cell presence.Furthermore, the low radscore group displayed significantly higher TNFRSF18 expression, while the high radscore group exhibited significantly elevated levels of HHLA2.TNFRSF18 is involved in immune activation, promoting T cell proliferation and enhancing anti-tumor immune responses, while HHLA2 plays a role in immune checkpoint regulation, suppressing T cell activation and contributing to immune evasion in the tumor microenvironment.Collectively, these comprehensive results provide further insights into the immune landscape associated with the radiomics signature.The observed activated immune microenvironment in the low radscore group and the immunosuppressive state in the high radscore group contribute to a better understanding of the underlying biological significance of the radiomics signature in predicting immunotherapy outcomes.
Our study has some limitations that need to be considered.Firstly, this is a retrospective study, which may inevitably lead to some selection bias in the collected information.Secondly, in the current study, we exclusively employed the median as the grouping threshold.Depending solely on a singular threshold, like the median, for patient stratification as a biomarker might prove insufficient to accommodate the intricacies of clinical practice [39].Thirdly, the methods for assessing MSI status differed between the cohort from Center 1 and the TCIA cohort.In addition, despite being a multicohort study, this research had a relatively small sample size and relied on publicly available databases.Future prospective multi-center studies with larger sample sizes are warranted to further validate the findings of this study.
In conclusion, our study developed a novel radiomics signature for predicting MSI-H expression in GC patients, guiding immunotherapy and predicting clinical outcomes.The radiomics signature also unveiled distinct immune profiles between low and high radscore groups, underscoring their clinical relevance.These findings emphasize the potential of radiomics analysis as a noninvasive tool for tumor characterization and personalized treatment selection in GC, warranting further validation in diverse populations and clinical settings to establish its clinical implications.

MSI
Microsatellite instability GC Gastric cancer TCIA The Cancer Imaging Archive LASSO The Least Absolute Shrinkage and Selection Operator

Fig. 2
Fig. 2 Radiomics feature selection by using the least absolute shrinkage and selection operator (LASSO) logistic regression.(a) The selection of tuning parameter (λ) in the LASSO model used 10-fold cross-validation via minimum criteria.The AUC curve was plotted versus log (λ).(b) LASSO coefficient profiles of the radiomics features.A vertical line was plotted at the optimal λ value, which resulted in 9 features with nonzero coefficients

Fig. 3 Fig. 4
Fig. 3 Receiver operating characteristic curves (ROC) for different models in the training (a), and validation cohorts (b); Calibration curves for the radiomics signature in the training (c), and validation cohorts (d); Radscore of different subtypes in the training (e), and validation cohorts (f).MSI-H = microsatellite instability-high, MSS = microsatellite stable

Fig. 5
Fig.5 The heatmap of the clinical-and immune-related molecular landscape.From the top to the end, there are five categories, encompassing clinical characteristics, immune cells, B7-CD28, TNF superfamily, and other immune-related molecular landscapes

Fig. 6
Fig. 6 Immune cell infiltration (a), proportion of microsatellite instability (MSI) status (b), expression status of TNFRSF18 (c) and expression status of HHLA2 (d) in different radscore groups

Table 1
Characteristics of patients used to evaluate MSI status Note: Except where indicated, data are number (%) of patients.MSI-H = microsatellite instability-high, MSS = microsatellite stability, CA 72 − 4 = carbohydrate antigen 72 − 4, CA 19 − 9 = carbohydrate antigen 19 − 9, CEA = carcinoembryonic antigen * Data in parentheses are interquartile range ROC curves for MSI-H positive expression are shown in Fig. 3a and b, and the results of the Delong test are presented in Supplementary TableE3.The analysis

Table 3
Multivariable logistic regression analysis of predictors for MSI status

Table 4
The specific performances of models for evaluate MSI status Note: Data in parentheses are 95% CI.MSI = microsatellite instability, AUC = area under the receiver operating characteristic curve, CI = confidence intervals